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High-performance, integrated spectrometers operating in the far-infrared and sub-millimeter promise 
to be powerful tools for the exploration of the epochs of reionization and initial galaxy formation. 
These devices, using high-efficiency superconducting transmission lines, can achieve the performance 
of a meter-scale grating spectrometer in an instrument implemented on a four-inch silicon wafer. Such 
a device, when combined with a cryogenic telescope in space, provides an enabling capability for 
studies of the early universe. Here, the optical design process for Micro-Spec (//-Spec) is presented, 
with particular attention given to its two-dimensional diffractive region, where the light of different 
wavelengths is focused on the different detectors. The method is based on the stigmatization and min- 
imization of the light path function in this bounded region, which results in an optimized geometrical 
configuration. A point design with an efficiency of ~ 90% has been developed for initial demonstra- 
tion, and can serve as the basis for future instruments. Design variations on this implementation are 
also discussed, which can lead to lower efficiencies due to diffractive losses in the multimode region. 

OCIS codes: (050.1940) Diffraction; (120.6200) Spectrometers and spectroscopic instrumen- 

tation; (220.2740) Geometric optical design; (300.6340) Spectroscopy, infrared. 


1. Introduction 

Far-infrared (IR) and sub-millimeter (15 /xm-1 mm) 
spectroscopy provides a powerful tool to probe a 
wide range of environments in the universe. In the 
past thirty years, many space-based observatories, 
such as Spitzer, Herschel, and the Infrared Space 
Observatory (ISO) [1-3], have opened the far-IR 
window to the universe, revealing rich line and con- 
tinuum spectra from objects ranging from inter- 
planetary dust particles to major galactic mergers 
and young galaxies in the early universe. Discov- 
eries made by these observatories have provided 
unique insights into physical processes leading to 
the evolution of the universe and its contents. How- 
ever, many questions remain unanswered regarding 
the very first moments of the universe, as well as 
galactic, stellar and planetary formation. 

Micro-Spec (/i- Spec) is proposed as a novel tech- 
nology concept to enable new discoveries in the far- 
IR spectral range, fi- Spec will be a high-sensitivity, 
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direct-detection spectrometer operating in the 450- 
1000-/tm regime. The light from the telescope is 
coupled to the spectrometer chip via a broadband 
slot antenna (Fig. 1) and directed to the power 
splitters and delay network realized in microstrip 
transmission lines. These microstrip elements are 
connected to a planar waveguide multimode region 
via transitions, which also serve as internal feed 
apertures. The planar waveguide multimode region 
has two internal antenna arrays, one for transmit- 
ting and one for receiving the spectrum as a func- 
tion of wavelength. Absorber structures lining the 
multimode region are used to control power emit- 
ted into large angles or reflected from the receiver 
antenna array. Finally, each receive antenna in the 
multimode region is coupled to a bank of order- 
sorting filters terminated in microwave kinetic in- 
ductance detectors (MKIDs) for readout. The in- 
strument, with a resolution 1Z ~ 1200, will be 
integrated on a four-inch-diameter (100 mm) sil- 
icon chip. This reduction in size is largely real- 
ized through the use of single-mode microstrip de- 
lay lines, which can be relatively compactly routed 
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Fig. 1: Layout of the fi - Spec module. The power is coupled into the instrument using a broadband 
antenna (left). It is then transmitted through a low-loss superconducting transmission line to a divider 
and a phase delay network, which creates a retardation across the input to the multimode region (in light 
blue). The feed horns will radiate a converging circular wave, which will concentrate the power along the 
focal surface, with different wavelengths at different locations. The outputs are connected to a bank of 
order-sorting filters to disentangle the various orders. (Color online) 


on the silicon wafer. Relative to counterparts re- 
alized in free space, the required transmission line 
lengths are more compact by a factor of the me- 
dia’s effective index. For all these reasons, /x-Spec 
can become an important instrument capability un- 
der the low background conditions realized in bal- 
loon and space-borne platforms such as the Space 
Infrared Telescope for Cosmology and Astrophysics 
(SPICA) [4]. The spectral range of the implementa- 
tion studied here is set by the materials employed 
for the transmission line and detectors. Applica- 
tion of this approach is presently confined to wave- 
lengths A > 250 /im by the gap frequency of avail- 
able low-loss superconductors. These include nio- 
bium (Nb) and niobium-titanium nitride (NbTiN) 
for the transmission line structures, and molybde- 
num nitride (MoN) for the detectors. 

Two distinct electromagnetic modes are used 
in the fi - Spec transmission spectrometer: mi- 

crostrip and parallel-plate waveguide [5] (see Fig. 4, 
right). Microstrip is a single-mode transmission line 
structure which is utilized in its dominant quasi- 
transverse electromagnetic mode (TEM) and in this 
context is analogous to single-mode fiber optic. It 
is used to realize the electrical delay line network. 
The parallel-plate waveguide has a TEM symme- 
try and the region serves as a multimode spatial 
power combiner. Its function can be compared 
to a grating spectrometer [6], in which a plane 
wave is reflected from the grating and the phase 
of each partial wave scattered from the rulings is 
a linear function of position across the grating. \i- 
Spec thus differs from similar technologies by the 
order of processing of the light in the spectrome- 


ter. For instance, in a Rowland spectrometer the 
required phase retardation is generated by reflec- 
tion from the grating grooves [6-10], whereas in 
Z-Spec, which is an example of planar Rowland 
grating architecture, propagation occurs in parallel- 
plate waveguides [11-14]. A last comparison can be 
made with one-dimensional bootlace lenses found 
in microwave practice [15-18], which /i-Spec builds 
upon for sub-millimeter wave applications. 

In this work the detailed optical design the /x- 
Spec transmission spectrometer is presented. We 
begin in Section 2 with an outline of the multimode 
power-combiner region optical design, then describe 
the internal two-dimensional feed structures used 
to couple from microstrip lines to this multimode 
region in Section 3, and finish by bringing these 
elements together to compute the modeled response 
of the entire antenna array in Section 4. 

2. Detailed spectrometer design 

The performance of the spectrometer is optimized 
by the selection of the x and y coordinates of each 
radiator and a differential phase shift on each de- 
lay line. In the case of the Rowland spectrometer, 
two points on the focal surface can be chosen, and 
wavelengths assigned to them [11]. The x and y po- 
sitions of the grooves in the grating can be selected 
to provide zero phase errors at each of these stig- 
matic points. A similar constructive approach can 
be used in the case of /i-Spec, where each emitter 
has an additional degree of freedom, the phase shift 
at the radiator, which allows generating a third stig- 
matic point. In the case of the Rowland grating, the 
selection of wavelengths at the two stigmatic points 
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determines the dispersion of the grating. 

In fi- Spec, the selection of the wavelengths of 
three stigmatic points determines both the disper- 
sion and the linearity of the dispersion of the sys- 
tem. A constraint on the selection is required 
if linear dispersion is desired. In the case pre- 
sented here, the stigmatic points and their asso- 
ciated wavelength were chosen to provide relatively 
linear dispersion and good imaging over the wave- 
length range of interest. The positions and wave- 
lengths of the stigmatic points are free parameters 
and can be selected to optimize aspects of the in- 
strument performance. A reasonable choice, which 
is adopted here, is to select points that minimize the 
overall root-mean-square (RMS) phase error for all 
spectrometer outputs (Fig. 3). 

In the prototype version of the instrument de- 
scribed here, in order to fit four spectrometers in a 
100-mm silicon wafer, each in a 25-mm x 25-mm re- 
gion, the maximum radius allowed is R — 1.25 cm. 
This choice is not optimal but enables several de- 
signs to be experimentally studied within a single 
mask set and the development of the required fab- 
rication process. The overall design was therefore 
carried out with the parameters shown in Table 1, 
where values for an initial design (A) and its opti- 
mized version (B) here adopted as the baseline are 
shown. 

With these parameters, the design procedure de- 
scribed was implemented in MATLAB through the 
goal attainment method of Gembicki [19]. This 
method was chosen because its implementation en- 
ables the resolution of vectorial equations subject to 
non-linear constraints which arise in the case pre- 
sented in this paper. The following set of equations 
describes the propagation of light through the spec- 
trometer, 


F m {x) = F^, m = l,...,3N e , (1) 

where F m represents the electrical delay which must 
be equal to i 7 ^, an integer multiple of the wave- 
lengths, at the three stigmatic points. N e repre- 
sents the number of radiators and 3 N e equations 
are needed to calculate the electrical delay of all 
radiators at each of the three stigmatic points. The 
vector of the unknowns, x, consists of the coordi- 
nates X{ and yi of the N e emitters’ centers, and the 
electrical path lengths in silicon, i??, for each feed 
horn delay line. Constraints are necessary to limit 
the search of the solution in the correct parame- 
ter space. In fact, each of the systems equations 
is of order 2 and the search of the solution was 
restricted within Ai/40 from the geometric optics 
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Fig. 2: Simplified representation of the grating 
geometry. On the left side three radiators can be 
seen, which point to the blaze point, B. The path 
from each radiator’s phase center to the first and 
second stigmatic points are also indicated by solid 
lines for the array’s i-th and center reference feed. 


solution for the grating solution (the 2 R dotted arc 
on the left side of Fig. 2). Using a numerical as 
opposed to an analytical solution thus allows us to 
analyze and incorporate the channel isolation and 
higher-order diffractive response into the optimiza- 
tion (see Sec. 3). Indeed, while the feed phase cen- 
ters reside on the position of the analytical solution, 
it will be shown that the array response is not uni- 
form in angle. 

All quantities in Eq. (1) are defined as follows: 


F(x) 


r>e I pe pe 

r>e I pe pe 

pe jl pe pe 

+ K Bi k b 


i = 1, ..., N e 


is the light path function matrix (3N e x 1), whose 
first 2 N e rows are associated with stigmatic points 
1 and 2 respectively, and the last N e rows with 
the blaze point B. All quantities in the F(x) ma- 
trix represent electrical lengths in the transmission 
medium. For example, the terms in the first N e 
rows associated with stigmatic point 1 are indicated 
below: 


• R\ is the electrical path length of the i-th sil- 
icon microstrip delay line; 
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Table 1: Spectrometer parameters for the initial (A) and optimized design (£>). Parameters without such 

a specific indication are invariant for either design. 


Minimum frequency 

h = c/Aoi = 450 GHz a 

Maximum frequency 

f 2 = c X o2 = 650 GHz 

Silicon, relative permittivity 

e r = 11.7 

Maximum wavelength in Si 

Ai = A 0 i / y/s r Hr — 195 /im 

Minimum wavelength in Si 

A 2 = A 02 / y/^rllr = 135 /iUl 

Blaze point wavelength in Si 

A b = VA 1 A 2 = 164 /im 

Antenna physical spacing, central array elements 

p(A) = 179 /im 


p(B) = 161 /im 

Design order of the grating 

M = 1 

Number of emitters 

N e = 65 

Number of receivers 

N r (A) =47 


N r (B ) =65 

Resolving power 

n = 6 5 

Multimode region, physical radius 

R(A) = 1.25 cm 


R{B) = 1.05 cm 


a Note: speed of light in vacuum, c = 3.00 x 10 8 m/s. 


•R e u = 


ne 


V ( R sin ]3i - Xip + ( R cos /?i - y;) 2 
x ^ /e r fi r , is the electrical path length of t 
i-th feed horn from stigmatic point 1 in the 
silicon parallel-plate waveguide multimode 
region; 


F* 


(i ~ M\ ol 
(i - ^tl) MX o2 
(i ~ Rf 1 ) MX oB 


: = l, ...,N e 


• R\ = 


(1 + sin /3i) 2 + cos 2 /?i 


X yj S r fl r , 

is the reference electrical path length of stig- 
matic point 1 from the central feed horn (i = 
(N e + l)/2), in the silicon multimode region. 


A 0.45-/im-thick monocrystalline silicon dielectric 
substrate is employed on the chip. The electrical 
length in the parallel-plate waveguide region is com- 
puted from the product of the index of refraction in 
the medium and the physical length defined by the 
geometry in Fig. 2. For example, R\ { = R\%^e r ii r ^ 
where e r and fi r « 1 are the effective relative per- 
mittivity and permeability specifying the propa- 
gation in the guiding medium. The parallel-plate 
waveguide conductor and ground plane are well ap- 
proximated and modeled as perfect electric conduc- 
tors (PEC). A surface reactance was assigned to 
the superconductor (niobium) to account for the 
kinetic inductance and model its influence on the 
impedance and propagation constant in the mi- 
crostrip transmission line structures. The super- 
conducting material was modeled as a 0.2-/im-thick 
layer with a London penetration depth, A^ = 95 
nm [20]. Similar definitions apply to stigmatic point 


where M — 1 represents the order of the grat- 
ing chosen for this design and A 0 the free-space 
wavelengths of the three stigmatic points. Here, 
the presence of subscript “o” signifies the radiation 
wavelength in free space and without indicates in 
the dielectric medium as defined in Table 1. In or- 
der to maximize the instrument efficiency, the emit- 
ters are tilted by an angle a so as to point to the 
center of the focal surface (Fig. 2): 


a.i — — arctan 


Vi A 
Xi\ + R ) 


( 2 ) 


A particular solution is found when placing the 
emitter antennas on the 2 R circle, using an axisym- 
metric configuration, and placing the blaze point on 
the optical axis. This provides a partial diagonal- 
ization of the problem, as in [15]. Under these hy- 
potheses, indeed, the last N e equations of (1), which 
refer to the blaze point, can be explicitly rewritten 
as shown below: 


m + 


y(i?sin/3 s - + {Rcos^b ~ yi^ ~ Rb 

N e + l 


X yj £ r fl r I i 


MX 


oBi 


X 


(3) 
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where (3 b is the angle between the y axis and the 
line from the center of the Rowland circle to the 
blaze point (Fig. 2). Here (3b = tt/ 2, which yields: 


Rf + 


= ^ 


(R - x,y- + yf - Rb 


X yj £ r fl r 


(4) 


Now, the terms y (R — xi) 2 + y 2 express the dis- 
tance of the z-th feed horn from the blaze point. 
Because the emitters lie on the 2 R circle, this quan- 
tity is exactly equal to 2 R. Also Rb is equal to 
2 R given the axisymmetric configuration of the in- 
strument and the particular location chosen for the 
blaze point. Therefore, these two terms cancel out 
and one readily obtains: 


R\ = M X oB 


N e + 1 


(5) 


the geometry of a Rowland spectrograph as a start- 
ing point, with the radiators on a surface of radius 
2 R, and the detectors on a circle of radius R at a 
distance of 2 R from the center of the radiators. Fi- 
nally, the spacing of the radiators is a free param- 
eter and will typically be chosen to limit diffrac- 
tion losses. We begin with a first radiator near the 
center of the radiator circle. Next, we choose an 
adjacent radiator and perturb its position and dif- 
ferential phase, which results in arriving at each of 
the three stigmatic points. This procedure is con- 
tinued for the radiators on each of the delay lines. 
This results in zero phase error at the stigmatic 
points by construction, low RMS phase errors be- 
tween the stigmatic points, and diffraction-limited 
performance in the multimode region over the de- 
sign spectral range. The RMS phase error, (-Prms > 
is illustrated in Fig. 3 for design A and it was cal- 
culated by assuming a uniform weight for each ra- 
diator at each point on the focal plane: 


which is a linear relation between the electric length 
in silicon, i??, and the feed horn number, z. 

To find the precise spacings of the emitters along 
the grating circle needed in order to satisfy the stig- 
matic equation for the pair of symmetric off-axis 
stigmatic points at specified high and low wave- 
lengths, our numerical solution method shows that 
the spacing between adjacent antennas is not a con- 
stant but grows quadratically as a function of posi- 
tion varying by ~ 13 ym from the nominal design 
value, p = 179 ym. This is the same behavior found 
by Rotman and Turner [15], who derived an analyt- 
ical solution for the analogous design of a bootlace 
lens, and found that the problem reduces to solv- 
ing an ordinary quadratic equation, from which one 
can determine the antenna positions. 

As mentioned in the beginning of this section, the 
diffractive solution can be found by restricting the 
algorithm to searching for radiator positions within 
a limited distance from the ideal Rowland circle. 
The constraint reads: 


\( (*< - -ft) 2 + yf - (2i2 + z) 


< Ai/40, (6) 


where z is the phase center correction provided by 
finite-element simulations (z = 15 /im, see Sec. 4) 
and Ai/40 ^ 5 ym is chosen to provide the upper 
and lower bounds. This is equivalent to a phase 
error of (2n/X 0 i) x (A o i/40) rad ~ 1/6 rad. 

Having chosen the nine parameters, i.e. , the co- 
ordinates of the stigmatic points and their wave- 
lengths, we proceed to determining the locations 
and phase shifts of the radiators. We have chosen 


VRMsifij) = 


f^Ve 

E 

_ 2= 1 


[}Pij (y(^j)))‘ 


nV2 


(7) 


Here, cpij is the relative phase of each transmit- 
ter, (ip(0j)) — 0 is the relative phase of the cen- 
tral transmitter (this is zero by construction as 
the central radiator is used as a reference, i.e., 
F (N e + 1)/2 = 0), N e is the number of transmitters, 
and Qj represents the angle corresponding to each of 
the points in which the focal plane was discretized. 
The phase between the z-th transmitter and the j- 
th point on the receiver focal plane is defined as 
(fij — 2iiyJe r (i r /\ 0 j x Fi and A 0 j is the free-space 
wavelength. 

3. Antenna feed response 

Conceptually, each microstrip feed structure is a 
two-dimensional analog of an adiabatic feed horn. 
To achieve the desired performance, the feed needs 
to match impedance, modal shape, and angular res- 
olution over the design frequency band. An adia- 
batic impedance taper [21] for the input microstrip 
line width to the output is achieved via a Hecken ta- 
per [22], which is a near-optimum matching section 
in the sense that it approaches an “equal-ripple” 
response without the use of abrupt discontinuities. 
These properties lower its sensitivity to fabrica- 
tion variations and tolerances. This feed section 
serves as a transition between microstrip lines and 
parallel-plate TEM waveguide regions. The latter 
is commonly referred to as the parallel-plate guide, 
and the microstrip planar waveguide model [5] can 
be used to compute its propagation and impedance 
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Fig. 3: RMS phase error between the average and the phase of all other transmit feed paths evaluated as 
a function of receiver position with the design parameters in Table 1 for design A. The RMS phase error 
vanishes at the three stigmatic points by construction and remains low in between. This ensures 
diffraction-limited performance over the power combiner’s design spectral range. The feed positions for 
the receive array in this design configuration are indicated by circular symbols in the figure. Addition of 
the radiators indicated by crosses can be used to improve the throughput in future implementations. 


properties. The transformer ratio, effective dielec- 
tric constant, and desired return loss set the neces- 
sary length of the impedance taper and were syn- 
thesized from readily available analytical forms and 
the response numerically simulated. A diagram of 
the feed aperture and spacing geometry is provided 
in the insert of Fig. 4 (left). In illuminating the 
multimode dielectric-filled parallel-plate waveguide 
region, the microstrip outputs play the role of the 
antenna apertures and need to be sized to minimize 
diffraction at the upper band edge. 

When used as a two-dimensional antenna ele- 
ment, the aperture geometry not only sets the ra- 
diation pattern but also plays a role in determin- 
ing the return loss [23]. This can be seen by the 
following consideration: the angular acceptance of 
the structure is constant over the band since the 
wavelength in the medium is less than the effec- 
tive aperture width, w. Considering the scaling 
of the reflection amplitude as a function of feed 
size [24], the desire for a low return loss at the 
low end of the band is in tension with the needs 
to limit diffraction (w < \ 0 / y/s r fir) and efficiently 
radiate (w > \ 0 /2 y Je r ii r ). This behavior effectively 
determines the smallest aperture size which can ef- 
fectively radiate [25] . The reflection from the end of 
the microstrip taper scales as a power of the wave- 
length over the microstrip width in two dimensions 
and is larger than that arising from the impedance 
taper, which can readily be made sub-dominant by 
increasing the transformer’s length. Thus, to ad- 
dress the reflection from the aperture, the output 
geometry of the aperture was apodized by adding 


a constant radius flare [26]. The resulting radiator 
element return loss, > 20 dB, is achieved over a 
2:1 bandwidth. Over the wavelength range where 
2 p < A, only the first-order diffraction occurs [27]. 
The influence of higher-order diffraction is miti- 
gated through the feed apodization, the coupling 
geometry of arrays, and the use of absorbing bound- 
aries on the multimode region walls. 

4. Antenna array response 

An efficient and accurate expression is desired for 
the evaluation of the structure’s response during 
numerical optimization. To evaluate the antenna 
array’s response, a classical phased array analysis 
approach is used. The emitters and receivers in 
the multimode region can be considered as a pair 
of opposing 2-dimensional antenna arrays. By con- 
sidering the Poynting vector as a function of spa- 
tial dimension, d, the following expression for the 
electric field is adopted for the computation of the 
array’s response: 


N e 


EM) = ]T 


Eq(0) 


2=1 


1+ (“) 


21 (d~l)/4 




(8) 


Here, r is the separation distance vector between 
the center point of the z-th emitter and a point 
in the multimode region, Eq is the E-field ampli- 
tude from the beam pattern angular distribution 
(9) at the center of the z-th feed horn (Fig. 4), 
w is the effective width of the feed horns, k = 
( 27i^/s r fi r /\ 0 ) • is the wavenumber vector in the 
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Parallel-Plate Waveguide 


0.2jL/m Nb 


0.2jLvm Nb 


Vacuum 


0.45jL/m Si 


Microstrip Line 


Fig. 4: (Left) The computed feed horn angular response at 430 GHz. The response of an emitting feed 
(aperture size 73 /im) is evaluated in the far field (||r|| ^ 6 A = 438 fim) and normalized to the magnitude 
of the E field at zero degrees. The feed array geometry is provided in the figure insert. The feeds phase 
center is indicated by a black filled circle. (Right) Cross section of the parallel-plate waveguide (upper) 
and microstrip (lower) transmission line geometries. The electric field is depicted with arrows to indicate 
the dominant modal symmetry. An adiabatic transition serves as a mode converter between these modes 
and defines the feed’s angular response in the multimode power combined region. 


medium, and is the unity vector normal to the 
z-th emitter. The phase c pij includes the contri- 
butions arising from the feed aperture’s phase, the 
phase delay, and the reference to the central emit- 
ter. 

The field amplitude’s dependence on radial sep- 
aration is physically motivated by the behavior of 
the fields in transitioning from the near to the far 
field [28] . Close in to the feed aperture the fields are 
essentially constant, and in going to large distance 
the amplitude decreases with distance. Physically 
useful limiting cases to consider are d — 3, which 
yields a scaling of the E- field inversely with distance 
for a spherical wave, and d — 1, which yields a scal- 
ing independent of distance for a lossless transmis- 
sion line. For d — 2, the case of primary interest 
here, it is important to note that the E-field am- 
plitude scales as the square root of the separation 
distance, ||r||, and is anticipated for a cylindrical 
geometry from the asymptotic form of the Hankel 
functions. This expression conserves power flow, is 
analytic over the parameter range of interest, and 
is a computationally efficient representation of the 
structure’s behavior during synthesis of the array 
response. The functional form and scaling indicated 
in Eq. (8) were validated by High-Frequency Struc- 
ture Simulation (HFSS) finite-element simulations 


of the structure’s response with distance. These 
simulations were used to derive the feed aperture’s 
effective width, re, and to determine the location of 
the phase center for radiation emitted by the struc- 
ture. 

The power is then computed from 

i / \ 1/2 

P(r, 0) = — E*(r, 0) ■ E(r, 0), r, = (^J Vo , 

(9) 

where r] 0 = 377 ft is the wave impedance of free 
space. Figure 5 illustrates the power distribution at 
450 GHz, as was calculated from Eqs. (9) and (8). 
On the left side, it can be seen that the power peaks 
at the first stigmatic point as required, and then 
it quickly drops to zero along the rest of the fo- 
cal plane. At separations from the feed elements 
greater than ~ 2 mm, the field transitions to the 
far field and the electric field amplitude is observed 
to asymptotically fall as ||r|| -1 ^ 2 . Simulations show 
the power is conserved throughout the multimode 
region. In this example, the ratio of the power emit- 
ted by the feed horns to the power received by the 
antennas is equal to ^ 30%. Most of the light at 
this wavelength is diffracted at angles that miss the 
receiver array, and in order to control reflections, 
are absorbed on walls on the two sides of the mul- 
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x (cm) x (cm) 

Fig. 5: Power distribution in the multimode region at a frequency of 450 GHz. On the left (Design A), 
the first-order peak is visible at an angle 6 « 10°, whereas the second-order diffraction peak shows up at 
6 ~ —66° (9 is defined in Fig. 2). In the optimized version (Design B) on the right, the higher-order 
diffraction peak has now disappeared and the throughput is as high as 87%. (Color online) 


timode region and between each receiver antenna. 
In addition, a second-order diffraction peak, visi- 
ble in lower left hand quadrant near the emitters, is 
responsible for part of this loss. 

Improved performance can be achieved with dif- 
ferent trade-offs between the antennas and the ar- 
ray geometry, which increase the throughput and 
reduce the effects of diffraction. If the two stig- 
matic points are moved away from the blaze point 
on the optical axis along the Rowland circle by 
4.2%, the instrument reaches its highest efficiency 
of 87% when the following applies: the radius is 
decreased by 16% ( R — 1.05 cm), the number of re- 
ceivers is increased by 18 ( N r = 65), and the emitter 
pitch, p, is set to 161 pm. The resulting increase in 
the receiver array solid angle is better matched to 
capture the transmitted beams given their apodiza- 
tion function. This is the primary reason for the 
multimode coupling-efficiency increase observed in 
this example. As an aside, this change in array ge- 
ometry reduces the magnitude of the second-order 
diffraction peak and moves it out of the range of 
angles of interest. For additional detail, see Fig. 5, 
rightmost inset. 

Fig. 6(a) shows the emitter antenna array’s nor- 
malized power response evaluated at the three stig- 
matic wavelengths for the optimized multimode re- 
gion design. Table 2 summarizes the coupling ef- 
ficiency between receiver and transmitter antenna 
arrays for the two design variations considered in 


this section. The inter-channel power isolation is 
presented in Fig. 6(b) for the center channel of 
the optimized design. The observed finite isolation, 
< 0.48, arises from the design goal to fully sample 
the spectrum and maximize the structure’s optical 
throughput. The resulting signal correlation can 
be addressed in calibration of the spectrometer re- 
sponse. Finally, the RMS phase error now reaches 
a maximum of 0.011 at /? = 75.8° and 104.3°, and 
vanishes at the three stigmatic points ((3 = 67.8°, 
90° and 112.2°). 

Table 2: Computed coupling efficiency (2p > A) 


Signal frequency Guide wavelength Coupling efficiency 
[GHz] [/xm] Design A Design B 


450 

Ai = 

= 195 

0.30 

0.91 

541 

A B z 

= 164 

0.30 

0.87 

650 

A 2 = 

= 135 

0.34 

0.55 


5. Conclusions 

A theory and design methodology has been de- 
scribed for p-Spec, a superconducting microstrip- 
based direct-detection spectrometer, in which the 
phase retardation is generated by propagation 
through transmission lines of different length. A 
proof-of-concept design was realized to demonstrate 
the feasibility of the instrument. The location of the 
three stigmatic points was chosen to minimize the 
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(a) (b) 


Fig. 6: Design B angular response, (a) The emitter antenna array’s angular power response is normalized 
by the total radiated power and plotted for the three stigmatic wavelengths. The directivity Z?o(A, 0) is 
the peak angular response of the array and is indicated on the plot, (b) Power isolation between the 
receiver channels computed for A 0 ~ 560 /im, where the power in each channel is normalized to the power 

in the central channel. 


overall RMS phase error for all spectrometer out- 
puts, whereas the location of both radiators and re- 
ceivers was found by setting the light path function 
equal to integer multiples of the three wavelengths 
associated with the three stigmatic points. The 
structure’s response was evaluated through an E- 
field model validated by finite-element simulations. 
The results show that, through a geometric opti- 
mization of the instrument, a maximum throughput 
of ^ 90% can be achieved in the multimode region. 
The attributes of this design approach, namely re- 
duced spectrometer envelope and high throughput, 
have the potential to enable a new class of space- 
born instrumentation which will reveal new details 
of the early universe. 
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